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摘 要 : 中 亚 地 区 常规 观测 站 点 稀少 , 需 借 助 星 载 高 光谱 AIRS 资料 分 析出 该 地 区 数值 预报 最 优 初 始 场 。 以 CRTM 


气 预报 业务 系统 中 的 适用 性 。 


中 输入 探 空 模拟 出 的 AIRS 辆 射 亮 温 为 参考 值 , 分 析 了 AIRS 观 测 亮 温 偏差 ,并 评估 了 ARS 卫星 资料 在 中 亚 数值 天 
结果 表明 :(1) 各 通道 模拟 所 选 站 点 上 空 


亮 温 最 大 正 偏差 的 平均 值 约 为 3.3 开 ,最 大 


负 偏 差 的 绝对 值 约 为 2.6K。(2) 多 个 站 点 平均 的 AIRS 观 测 辐射 亮 温 整体 咯 高 于 模拟 亮 温 ,其 概率 密度 分 布 比 单个 


站 点 更 加 接近 正 态 分 布 曲线 。(3) AIRS 的 同化 改善 了 RMAPS-CA 对 位 势 高 度 .温度 . 比 湿 等 高 空 要 素 的 预报 效果 ， 


并 未 改善 高 空 风 速 的 预报 。 对 各 个 要 素 ,AIRS 的 同化 改善 幅度 在 低层 较 高 层 大 。 同 化 后 ,位 势 高 度 温度 . 比 湿 和 
风速 的 预报 RMSE 分 别 小 于 20 gpm 2 K 8x10" kg*kg"' 以 及 5 meso 
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中 亚 地 区 位 于 欧 亚 大 陆 腹 地 ,是 北半球 最 大 的 
干旱 . 半 和 干旱 地 区 。 中 亚 地 区 泛 指 亚洲 中 部 ,及 伊 
朗 、 阿 富 汗 .巴基斯坦 .印度 的 部 分 地 区 也 应 属于 中 
亚 范 畴 '"。 受 中 纬度 西风 带 的 控制 ,中 亚 地 区 为 典 
型 的 大 陆 性 气候 ” ,冰雪 融 水 为 该 地 区 的 重要 水 源 
补给 ” ,研究 及 预报 中 亚 地 区 的 降水 及 温度 变化 等 
已 成 为 广大 气象 工作 者 的 研究 热点 中 。 然 而 ,由 于 
常规 观测 稀 玻 ,要 分 析出 该 地 区 数值 天 气 预 报 系 
统 的 最 优 初 始 场 , 仅 依靠 少量 的 常规 观测 资料 远 远 
不 够 ,还 需 借 助 高 分 辩 率 的 星 载 大 气 红外 探测 器 
AIRS 资 料 ”"。 

先进 的 红外 光谱 新 型 探测 设备 AIRS 是 美国 宇 
航 局 AQUA 卫星 上 携带 的 众多 仪器 之 一 ,是 能 采集 
2378 个 停 贷 的 热 红 外 辐射 的 高 光谱 红外 探测 仪 。 
采用 红外 光栅 分 光 技 术 ,高 光谱 分 辩 率 O/A = 
1200 ,每 天 全 球 覆 善 2 次 ,可 以 连续 在 3.7~15.4 um 
的 光谱 范围 内 观测 地 球 大 气 及 变化 。AIRS 的 扫描 
宽度 为 1650 km. AIRS 的 星 下 点 空间 分 辩 率 为 
13.5 km ,垂直 分 辨 率 1 km, 4$ 2.667 s 横 向 扫描 1 次 ， 
扫描 角度 +49.5° ,每 个 观测 视 场 包含 2378 个 光谱 取 
样 ,90 个 视 场 组 成 一 个 扫描 线 ,扫描 线 长 1650 km, 
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每 6 min 的 135 条 扫描 线 组 成 一 个 景 ,每 天 全 球 有 
240 个 景 的 观测 ,其 中 约 120 个 景 在 白天 ,而 另外 
120 个 景 在 夜间 。AIRS 的 大 多 数 信道 在 光谱 响应 和 
灵敏 度 方面 表现 出 优异 的 性 能 "。 但 观测 与 模拟 
辐射 之 间 的 系统 偏差 客观 存在 ， 使 得 其 推 广 应 用 仍 
然 面临 许多 问题 。 在 AIRS 直接 同化 应 用 到 数值 预 
报 业 务 系统 "之前, 须 将 卫星 观测 亮 温 进行 质量 
控制 ”和 偏差 订正 一 5 ,而 偏差 分 析 “是 偏差 订正 
的 第 一 步 。 国 内 分 析 AIRS 反 演 温 湿度 廓 线 偏差 的 
人 研究 较 多 '" ,分析 AIRS 探测 亮 温 偏差 的 研究 并 不 
ZL", An, AIRS 辆 射 亮 温 业 务 已 经 应 用 于 我 
国 数值 预报 中 心 的 数值 天 气 预 报 业 务 系统 ,并 提高 
了 数值 预报 准确 率 ”” ,但 在 中 亚 数 值 天 气 预报 业 
务 系 统 RMAPS-CA(Regional Meso-scale Analysis and 
Prediction System in Center Aisa) ji AFF JR AIRS 辆 射 
亮 温 资 料 的 同化 。 

鉴于 此 ,本 研究 首先 通过 与 CRTM (Community 
Radiative Transfer Model ) 2 模拟 亮 温 相 比 ,分 析 了 
AIRS 卫星 亮 温 的 偏差 分 布 特征 ,然后 用 美国 国家 大 
气 研究 中 ， 忆 (NCAR) 开 发 的 中 尺度 模式 WRF(ARW) 
V4.0 及 其 三 维 变 分 同化 系统 WRF-3DVAR™” ,以 
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2016 年 7 月 1 日 晴空 天 气 为 例 , 评 估 了 AIRS 卫星 资 
料 RMAPAS-CA 中 的 同化 效果 ,为 AIRS2 在 RMAPA 
S-CA 中 的 业务 同化 应 用 奠定 技术 基础 。 


1 资料 和 方法 


1.1 资料 

本 文 所 用 的 辐射 亮 温 资 料 是 AIRS 的 一 级 产品 ， 
可 从 GDAS(global data assimilation system ) 免 费 获取 
(ftp://ftpprd.ncep.noaa.gov/pub/data/necf/com/efs/prod/) o 
AIRS 每 天 有 00UTC ,06UTC, 12UTC,18UTC 共 4 个 时 
次 的 二 进 制 BUFR 格式 数据 文件 ,每 个 时 次 的 数据 
为 将 每 个 时 次 前 后 3 的 数据 拼接 到 一 起 生成 的 文 
件 。 在 日 常 的 数值 天 气 预报 业务 同化 中 ,常用 到 的 
AIRS 的 281 个 辐射 亮 温 通 道 , 故 在 本 次 偏差 分 析 及 
适用 性 评估 中 , 仅 对 其 中 的 281 个 通道 进行 研究 。 
本 文 所 用 的 探 空 资料 观测 时 间 为 2016 年 7 月 1 个 
月 ,每 日 观测 2 次 ,分 别 为 07:15 .19:15( 北 京 时 )。 本 
研究 用 于 驱动 区 域 数 值 预报 模式 的 全 球 背 景 场 为 
0.5° 的 GFS(Global Forecast System ) 。 
1.2 方法 

文中 用 到 的 辐射 传输 模式 CRTM (Community 
Radiative Transfer Model) 是 由 卫星 资料 同化 联合 中 
心 (Joint Center for Satellite Data Assimilation, JCS- 
DA ) 开 发 的 快速 辐射 传输 模式 , 它 是 针对 卫星 资料 
同化 而 开发 的 ,对 云雨 条 件 下 的 卫星 观测 的 模拟 计 
算 能 力 较 强 。CRTM 模式 主要 由 前 向 模块 、 切 线性 
模块 伴随 模块 以 及 K 和 矩阵 模块 等 4 个 模块 组 成 , 主 
要 功能 模型 包括 下 垫 面 发 射 反射 模型 `. 气 溶胶 吸收 
反射 模型 \ 云 吸收 散射 模型 .气体 吸收 模型 和 辐射 
传输 方程 求解 模块 。 本 研究 利用 WRF 模 式 对 研究 区 
域 和 时 间 进 行 预报 ,得 到 CRTM 模 拟 亮 温 所 需 的 地 表 
类 型 和 云 量 等 相关 参数 ,然后 参考 CRTM 中 
LOAD_ATM.in 文 件 中 的 温度 .气压 等 7 项 要 素 的 格 
式 将 TTAA 格 式 的 探 空 报 文 资料 处 理 成 98 层 的 廓 线 ， 
进一步 输入 到 CRTM 模 式 中 ,模拟 出 AIRS 红 外 亮 温 ， 
并 以 此 为 真 值 ,分 析 通 过 与 AIRS 观 测 亮 温 的 偏差 。 

本 文 将 RMAPS-CA 的 母 网 格 区 域 作为 研究 区 
域 。 该 系统 使 用 的 中 下 度 模式 是 WRF(ARW)V4.0， 
模拟 采用 单 层 网 格 方案 (图 1)。 具 体 的 数值 模拟 方 
案 ” 见 表 1。 

由 于 GFS 中 首先 已 同化 了 AIRS 观 测 , 故 在 设计 
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图 1 研究 区 示意 图 


Fig. 1 Target region 


表 1 数值 模拟 方案 


Tab. 1 Numerical simulation schemes 


参数 名 称 参数 设置 
初始 条 件 和 边界 条 件 GFS (0.5°x0.5°) 
中 心 点 (43.55°N ,87.85°E) 
格 点 712 kmx532 km 
格 距 9 km 
垂直 层 数 50 
模式 层 顶 气压 10 hPa 
时 间 步 长 10s 
微 物理 过 程 方案 Thompson 
长 波 辐射 方案 RRTMG scheme 
短波 辐射 方案 RRTMG scheme 
积 云 对 流 方 案 Kain-Fritsch 
陆 面 过 程 方案 Noah Land Surface Model 
边界 层 方案 ACM2 


数值 试验 时 ,用 了 提前 6h 的 CFS 预报 场 作为 背景 
场 , 即 用 2016 年 6 月 30 日 18:00 的 CFS 预报 场 作为 
2016 年 7 月 1 日 00:00 的 背景 场 。NR 试验 中 不 同化 
任何 资料 ,而 AIRS 试验 同化 2016 年 7 月 1 日 00:00 
的 AIRS 资 料 , 依 此 类 推 。 


2 偏差 分 析 


分 析 中 亚 地 区 的 AIRS 亮 温 偏 差 ,首先 要 知道 该 
区 域 AIRS 资料 的 空间 分 布 情况 。 对 2016 年 7 月 1 
日 的 00UTC、06UTC、12UTC 和 18UTC4 个 时 次 的 
AIRS 辐 射 亮 温 观测 (图 2 中 的 X) 按 照 气压 (P) 属 性 
进行 分 层 , 各 层 对 应 气压 范围 见 表 2。 
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R2 对 AIRS 辐射 率 观测 的 气压 分 层 设 定 


国境 内 的 亮 温 观测 多 在 气压 大 于 700 hPa 以 下 的 低 


Tab. 2 Pressure layers setup of observed AIRS radiance 空 。 而 低空 亮 温 观测 由 于 受到 地 表 辐 射 .地 形 等 因 
分 层 序号 最 大 压 /hPa 最 小 气压 /hPa 素 的 影响 ,使 得 其 与 CRTM 模拟 的 探 空 亮 温 可 能 有 
i. beige 较 大 偏差 。 在 12UTC+3h 时 间 窗 口内 (图 2b)， 
p pa AIRS_AQUA 绝 大 多 数 观测 集中 在 中 亚 偏 西 地 区 。 
4 200 150 在 中 国 的 偏 西 偏 北 地 区 分 布 较 00UTC+3 h 时 间 窗 口 
5 300 200 内 的 AIRS 辐射 亮 温 观测 数量 更 稀 玻 ,也 主要 以 P< 
6 500 300 50 hPa 的 高 层 亮 温 观 测 为 主 (图 2c)。 在 18UTC+3h 
a si 时 间 窗 口内 ,AIRS_AQUA 扫描 带 有 2 个 , 偏 西 的 1 
8 1000( 模 式 层 底 ) 700 


个 扫描 之 扫 过 中 国 新 疆 、 甘 肃 和 青藏 高 原 大 部 ,而 


从 图 2 可 以 看 到 ,在 00UTC+3h 时 间 窗 口内 , 绝 ” 也 较 06UTC+3h 时 间 徐 口内 的 东部 扫描 带 整 体 向 东 
大 多 数 观测 集中 在 中 亚 偏 西 地 区 。 在 中 国 的 偏 西 。 偏 移 了 约 5 个 经 度 ( 图 24)。 


俩 北 地 区 分 布 稀 玻 , 且 主 要 以 P<50 hPa 的 高 层 亮 
温 观 测 为 主 ( 图 2a)。 在 06U 


鉴于 AIRS 辐 射 亮 温 观 测 在 中 亚 地 区 的 分 布 情 


TC+3h 时 间 窗 口内 ，” 况 (图 2), 结 合 RMAPS-CA 的 母 网 格 覆 盖 范 围 ( 图 


AIRS_AQUA 扫描 带 有 2 个 , 偏 西 的 1 个 扫描 带 扫 过 ”了 ,本文 选择 了 06UTC+3h 和 18UTC+3 h 2 个 时 间 
中 国 新 疆 和 西藏 的 大 部 分 区 域 ,而 偏 东 的 扫描 带 扫 ”窗口 内 窗 盖 中 国 西部 地 区 的 扫描 带 里 的 探 空 站 作 
过 内 蒙古 和 甘肃 直至 云南 的 中 国 中 部 一 带 , 且 在 中 为 代表 站 点 ( 表 1)。 与 此 同时 , 因 00UTC+3h 和 


(b) 2016-07-01 06 UTC 
sce ai 


wi: 


Bech 


Yast 
wt 


70°E 
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e 200 hPa<x<300 hPa 


80°E 90° 
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@ 50 hPa<x<100 hPa e 100 hPa<x<150 hPa @ 150 hPa<x<200 hPa 


图 


e 300 hPa<sx<500 hPa = @ 500hPa<x<700hPa ® x>700 hPa 


2 AIRS 辐射 亮 温 空间 分 布 


Fig. 2 Spatial distribution of AIRS brightness temperature (BT) 


202103.00066v1 


chinaXiv 


马 玉 芬 等 :AIRS FAT Se TEP DH XA 2E PT ee AE 


12UTC+3 h 2 个 时 间 窗 口内 AIRS_AQUA 扫描 带 的 
星 下 点 与 RMAPS-CA 的 母 网 格 中 心 偏离 较 远 ,可 能 
导致 模拟 的 探 空 亮 温 与 观测 亮 温 出 现 较 大 偏差 ,而 
这 种 大 偏差 的 亮 温 观测 在 业务 同化 时 也 会 被 WRF- 
DA 自动 剔除 , 故 不 选择 这 2 个 时 段 内 的 卫星 扫描 带 
内 的 探 室 站 点 进行 研究 。 

在 筛选 出 的 站 点 ( 表 3) 中 ,我 们 选择 在 站 点 经 
纬度 +0.2* 范 围 内 ` 离 观测 站 点 距离 最 短 .观测 时 间 
间隔 最 小 的 AIRS 辐射 亮 温 与 CRTM 模拟 的 探 空 亮 
温 对 比 进 而 分 析 AIRS 辆 射 亮 瘟 。 此 处 分 析 的 AIRS 
辐射 亮 温 , 仅 限于 在 数值 天 气 预报 业务 系统 同化 模 
块 中 常用 的 281 个 通道 。 受 限于 上 述 空间 配对 方 
法 ,筛选 出 了 表 2 中 的 12 个 站 点 ,分 析 了 2016 年 7 
月 1 一 31 日 ( 塔 中 站 加 密 观 测 实验 持续 时 间 ) 期 间 各 
个 站 点 的 整 月 辐射 亮 温 平均 偏差 (图 3 )。 

从 代表 探 空 站 点 的 AIRS 的 亮 温 偏差 (图 4) 可 
以 看 到 ,各 个 站 点 处 模拟 出 的 探 空 亮 温 数量 各 不 相 
同 , 均 小 于 预期 模拟 出 的 281 个 通道 数 , 模 拟 出 的 通 
道 数 量 最 大 值 191 在 伊 宁 站 ,最 小 值 65 在 塔 中 站 ， 
仅 塔 中 站 和 民 丰 站 模拟 控 空 亮 温 通道 数 在 100 以 
下 ,其余 10 个 站 的 通道 数 在 163~190 之 间 , 约 占 了 
所 有 281 个 亮 温 通道 数 的 58%~68%。 与 此 同时 ,各 
通道 模拟 探 空 亮 温 偏差 有 正 有 人 负 ,整体 偏差 范围 
在 -4.0K~4.0K。 各 站 的 最 大 和 最 小 偏差 不 同 , 且 最 
小 偏差 的 绝对 值 小 于 最 大 偏差 的 绝对 值 ,最 大 偏差 
的 平均 值 约 为 3.3, 最 小 偏差 的 绝对 值 约 为 2.6。 其 
中 , 明 丰 站 的 最 大 偏差 比 其 他 11 个 探 空 站 的 最 大 偏 
差 小 , 仅 为 0.6, 最 小 偏差 的 绝对 值 也 比 其 他 11 个 探 
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空 站 小 ,最 小 偏差 仅 为 -0.6( 表 1)。 

资料 偏差 符合 正 态 分 布 的 特点 是 三 维 变 分 资 
料 同化 的 必要 条 件 之 一 。 为 进一步 探讨 探 空 模拟 
亮 温 与 AIRS 观测 亮 温 的 偏差 ,我 们 绘制 了 亮 温 偏差 
的 概率 密度 分 布 直方 图 ,并 拟 合 了 其 高 斯 概率 密度 
分 布 曲 线 。GaussAmp 函数 的 一 般 形式 为 : 


f(x) =o + ae 8" 
式 中 :a.5 与 c 为 实数 常数 , 且 a>0。a 为 高 斯 概率 
密度 曲线 的 峰值 。! 为 峰值 对 应 的 横 坐 标 , 即 均 
数 。c 为 标准 差 。yo 为 概率 密度 曲线 的 基数 。 假 设 
一 般 的 正 态 总 体 在 某 个 特定 区 间 的 概率 总 体 用 
COD 来 表示 。 

从 代表 站 点 的 AIRS 的 亮 温 偏差 概率 密度 曲线 
基态 值 m 来 看 ,12 个 探 空 站 中 仅 伊 宁 、 库 车 喀什、 
RE 塔 中 站 的 模拟 AIRS 辐射 亮 温 偏差 概率 密度 曲 
线 的 基态 在 2.0 以 下 ,比较 接近 正 态 曲线 的 m=0。 另 
外 的 7 个 探 空 站 模拟 AIRS 辐射 亮 温 偏差 概率 密度 
曲线 与 正 态 分 布 相差 较 远 ,y 在 2.0 以 上 ,最 大 达到 
了 3.46, 在 阿勒泰 站 。 

由 图 4 可 知 ,从 峰值 a 来 看 ,12 个 探 空 站 的 模拟 
AIRS 辐射 亮 温 偏差 概率 密度 曲线 的 峰值 最 大 值 为 
24.77 ,在 阿勒泰 站 ;最 小 值 为 8.10 ,在 库 车 站 。 另 外 
的 10 个 探 室 站 中 , 库 车 、 民 丰 、 塔 中 站 的 峰值 也 较 
小 ,在 10 以 下 ,其 余 7 个 站 的 峰值 在 10~20 之 间 。 
从 均 数 2 来 看 ,12 个 探 空 站 的 模拟 AIRS FSS eH id 
差 概率 密度 曲线 的 均 数 在 -0.6~ 0.88 之 间 ,绝对 值 
均 小 于 1。 其 中 绝对 值 的 最 大 值 0.88 在 喀什 站 ,说 
明 喀 什 站 的 的 概率 密度 曲线 的 对 称 性 较 差 ,最 小 


表 3 代表 站 点 亮 温 值 及 相关 参数 


Tab.3 Mean bias of BT over the selected sounding sites and other relevant parameters 


序号 探 空 站 点 探 空 站 点 经 纬度 通道 数 最 大 偏差 最 小 偏差 基态 峰值 ” 均 数 ” 标准 差 
1 阿勒泰 (88.08°E,47.73°N) 179 3.8 -2.8 3.46 22.05 -0.35 0.13 
2 AIR (83.00°E ,46.73°N ) 179 4.0 -2.8 2.68 11.38 -0.28 0.35 
3 克拉 玛 依 (84.86°E ,45.63°N ) 179 4.0 -2.8 3.22 10.89 -0.29 0.34 
4 FF (81.33°E,43.95°N ) 190 3.0 -2.6 0.21 11.07 0.61 1.30 
5 乌鲁木齐 (87.65°E ,43.78°N ) 163 3.4 -3.4 2.01 13.00 -0.60 0.39 
6 阿克苏 (80.23°E,41.12°N) 186 3.2 -2.2 2.82 19.87 0.50 0.18 
7 库 车 (83.07°E ,41.72°N ) 171 3.4 -3.2 -0.03 8.10 0.44 1.72 
8 库尔勒 (86.13°E,41.75°N) 170 4.0 -2.8 2.69 15.19 0.43 0.24 
9 喀什 (75.75°E ,39.48°N ) 187 2.8 -2.2 1.89 13.19 0.88 0.59 
10 RE (82.72°E ,37.07°N) 77 0.6 -0.8 1.47 11.58 1.31 0.24 
11 和 | (79.93°E ,37.13°N ) 179 3.0 -2.8 1.07 9.92 -0.45 0.86 
12 塔 中 (83.63°E ,39.04°N ) 65 3.6 -3.0 0.84 8.23 0.75 1.15 
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图 3 代表 探 空 站 点 的 AIRS 的 亮 温 偏差 


Fig.3 BIAS of AIRS BT over the selected sounding stations 


值 -0.28 在 塔 城 站 和 克拉 玛 依 站 ,说 明 这 2 个 站 的 概 
率 密度 曲线 的 对 称 性 较 好 。 从 标准 差 c 来 看 ,模拟 
AIRS 辐 射 亮 温 偏差 概率 密度 曲线 的 标准 差 在 0.03~ 
0.19 之 间 。 其 中 , 民 丰 站 和 阿克苏 站 的 标准 差 最 
小 。 整 体 来 看 , 仇 宁 站 、 库 车 站 、 和 田 站 、 塔 中 站 模 
拟 AIRS 辐 射 亮 温 偏差 概率 密度 较 接近 正 态 分 布 , 阿 
勒 泰 .阿克苏 .库尔勒 站 模拟 AIRS 辐射 亮 温 偏差 概 
率 密度 与 正 态 分 布 相差 最 远 ,其 余 站 居中 。 

与 单个 站 点 的 偏差 概率 密度 分 布 相 比 ,12 个 探 
空 站 点 的 平均 偏差 概率 密度 分 布 更 加 接近 正 态 分 
布 曲线 ,其 概率 密度 曲线 的 基态 值 y=3.44, 峰 值 a~ 
100 ,均值 上 约 等 于 0.33 ,表示 多 个 站 点 平均 的 AIRS 


辐射 亮 温 观测 整体 略 高 于 CRTM 模拟 的 探 空 亮 温 
月 平均 正 偏差 为 0.33 人 。 标 准 差 c=1.74( 图 5)。 


3 适用 性 评估 


从 偏差 分 析 的 结果 可 以 看 到 ,虽然 AIRS 辆 射 亮 
温 的 平均 偏差 多 在 2K 以 内 且 偏 差 概率 密度 接近 正 
态 分 布 , 但 个 别 通道 的 月 平均 偏差 值 仍然 较 大 。 为 
一 步 分 析 AIRS 辐射 亮 温 资料 在 中 亚 区 域 的 适用 
性 ,我们 设计 了 一 组 AIRS 辐射 率 同化 的 敏感 性 实 
验 , 在 不 进行 资料 同化 的 控制 试验 (NR) 基 础 上 , 设 
计 了 AIRS 红 外 辐射 亮 温 资料 同化 试验 (AIRS ) 。 
鉴于 AIRS 辐射 亮 温 观测 的 空间 分 布 ,我 们 在 进 
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图 4 代表 探 空 站 点 模拟 的 AIRS 亮 温 偏 差 概率 密度 分 布 
Fig.4 Frequency density of AIRS BT over the selected sounding stations 


行 同化 效果 评估 时 , 仅 对 位 势 高 度 .温度 和 比 湿 等 
高 空 要 素 的 模拟 值 相对 于 模拟 区 域 的 探 空 值 的 偏 
差 (BIAS) 和 均 方 根 误差 (RMSE ) 进 行 定量 分 析 。 由 
于 假设 为 真 值 的 探 空 资 料 时 间 间 隔 为 12 hn, 故 我们 
对 模式 积分 初始 时 刻 ( 以 下 简称 00UTC) 及 积分 12h 
(以 下 简称 12UTC) 和 24 h( 以 下 简称 24UTC) 这 3 个 
时 次 的 高 空 要 素 BIAS 和 RMSE 分 别 进行 对 比分 析 ， 
比较 有 无 AIRS 辐射 率 同化 对 各 高 空 要 素 在 各 个 时 


次 预报 准确 率 的 影响 。 

(1) 从 位 势 高 度 场 预 报 BIAS (图 6a) 和 RMSE 
(图 7a) 来 看 。 在 00UTC ,NR 实验 中 ,位 势 高 度 预报 
BIAS 整 体 为 正 , 且 在 850 hPa 的 低层 达到 了 最 大 值 
约 120 gpm。 而 在 AIRS 实 验 中 ,200 hPa 以 下 的 低空 
位 势 高 度 预报 RMSE BE NR 实验 有 较 明 显 的 减 小 , 低 
空 的 位 势 高 度 正 BLAS 调整 为 负 , BLAS 绝对 值 减 
小 。 在 200 hPa 以 上 的 高 空 ,AIRS 同化 对 位 势 高 度 
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140'[ e BIAS 整 体 为 负 且 随 高 度 变化 不 大 。BIAS 绝 对 值 均 
a=100.26137 a 
120 + 均值 p=0.33166 在 20 gpm 以 内 , 仅 在 100 hPa 附近 的 高 空 为 正 BIAS， 


标准 差 c=1.7419 


0 
平均 偏差 用 


图 5 模拟 的 AIRS 亮 温 站 点 平均 偏差 概率 密度 分 布 
Fig. 5 Frequency density distribution of averaged BIAS of 
the simulated AIRS BT 


的 预报 为 负 效 果 。 在 12UTC,NR 实验 中 ,位 势 高 度 
预报 BLAS 整体 为 负 且 随 高 度 变 化 不 大 , 仅 在 200 
hPa 以 上 的 高 空 为 正 。BIAS 绝对 值 和 均 方 根 均 在 
20 gpm ANo 与 NR 实验 相 比 ,AIRS 实 验 500 hPa 以 
下 的 低空 位 势 高 度 预报 RMSE 有 所 减 小 ,其 上 各 层 
效果 不 明显 。 在 24UTC,NR 实验 中 ,位 势 高 度 预报 


RMSE 在 25 gpm 以 内 , 且 在 200 hPa 附近 RMSE 最 
小 。 在 AIRS 实 验 中 ,位 势 高 度 预 报 BIAS 整 体 也 为 
负 且 随 高 度 变化 不 大 。BIAS 绝 对 值 均 在 20 gpm 以 
内 , 仅 在 100 hPa 附近 的 高 空 为 正 BIAS,RMSE 在 25 
gpm 以 内 , 且 在 200 hPa 附近 RMSE 最 小 。 与 NR 实 
验 相 比 ,AIRS 同化 对 24UTC 的 位 势 高 度 的 预报 整体 
为 负 效 果 。 综 合 来 看 ,对 位 势 高 度 的 预报 ,AIRS 同 
化 在 00UTC 和 12 UTC 时 次 为 正 效果 ,在 24UTC 为 
负 效 果 , 同 化 效果 在 各 个 时 次 在 500 hPa 以 下 的 低 
层 整体 较 优 。 

(2) 从 高 空 温度 场 预报 BIAS( 图 6b) 和 RMSE 
(图 7b) 来 看 ,在 00UTC,AIRS 实 验 各 层 温 度 预报 BI- 
AS 绝对 值 和 RMSE 均 比 NR 实验 小 ,而 在 12UTC 和 
24UTC 整 体 较 NR 大 。 从 BIAS 来 看 ,在 各 个 预报 时 
Uk, AIRS 辐射 率 同 化 使 得 300 hPa 以 上 的 高 空 各 层 
温度 整体 变 暖 ,其 下 变 冷 ,调整 幅度 最 大 在 850 hPa 
以 下 的 低层 。 同 化 AIRS 后 ,各 层 温 度 偏 差 绝 对 值 和 
RMSE 均 在 2K 以 内 。 
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图 6 高 空 要 素 预 报 侦 差 


Fig. 6 Forecast bias of spacial meteorological variables 
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图 7 高 空 要 素 预 报 均 方 根 误差 


Fig. 7 Forecast RMSE of spacial meteorological variables 


(3) 分 析 比 湿 预 报 BIAS( 图 6c) 和 RMSE ( Al 
7c)。 比 湿 (specific humidity) 是 指 在 一 团 湿 空气 中 , 
水 汽 的 质量 与 该 团 空气 总 质量 (水 汽 质 量 加 上 干 空 
气质 量 ) 的 比值 。AIRS 辐射 率 同 化 在 各 个 时 次 对 比 
湿 场 预报 都 产生 了 不 同 程度 的 改进 。 由 于 NR 实验 
和 AIRS 实 验 中 ,明显 的 BIAS 均 出 现在 300 hPa 以 下 
的 低空 , 故 以 下 主要 分 析 300 hPa 以 下 比 湿 的 预报 
BIAS 和 RMSE。 与 NR 实 验 相 比 ,AIRS 实验 在 各 个 
预报 时 效 对 比 湿 场 预报 的 RMSE 整体 较 小 ,说 明 
AIRS 同化 对 比 湿 预 报 整体 为 正 效 果 。 同 化 后 , 比 湿 
场 的 RMSE 在 各 层 均 小 于 8x10”kg'kg 。 

(4) 分 析 高 空 风速 预报 的 BIAS( 图 6d) 和 RMSE 
(图 7d)。 同 样 的 ,在 各 个 预报 时 次 ,AIRS 的 同化 对 
850 hPa 以 下 的 低层 的 风速 调整 幅度 最 大 ,但 同化 
并 未 产生 明显 的 正 效果 。 


4 ”结论 与 讨论 


文章 分 析 了 AIRS 观测 卫星 亮 温 的 偏差 分 布 特 
征 ,评估 了 AIRS 卫星 资料 同化 效果 ,得 到 如 下 结论 : 


(1) 在 所 选 的 12 个 探 空 站 上 空 ,各 通道 模拟 探 
空 亮 温 最 大 偏差 的 平均 值 约 为 3.3 ,最 小 偏差 的 绝 
对 值 约 为 2.6K。 伊 守 站 、 库 车 站 、 和 田 站 、 塔 中 站 模 
拟 AIRS 辐射 亮 温 偏差 概率 密度 较 接近 正 态 分 布 , 阿 
WAS 阿克苏、 库尔勒 站 模拟 AIRS 辐 射 亮 温 偏差 概 
率 密度 与 正 态 分 布 相差 最 远 ,其 余 站 居中 。 

(2) 与 单个 站 点 的 偏差 概率 密度 分 布 相 比 ,12 
个 探 空 站 点 的 平均 偏差 概率 密度 分 布 更 加 接近 正 
态 分 布 曲 线 ,其 概率 密度 曲线 的 基态 值 y=3.44 , 峰 
值 o=100 ,均值 ! 约 等 于 0.33 ,表示 多 个 站 点 平均 的 
AIRS 辐射 亮 温 观测 整体 略 高 于 CRTM 模拟 的 探 空 
亮 温 ,月 平均 正 偏差 为 0.33 K。 标 准 差 c~1.74。 

(3) AIRS 的 同化 改善 了 RMAPS-CA 对 位 势 高 
BE .温度 . 比 湿 等 高 空 要 素 的 预报 效果 ,并 未 改善 高 
空 风 速 的 预报 。 对 各 个 要 素 ,AIRS 的 同化 改善 幅度 
在 低层 较 高 层 大 。 同 化 后 ,位 势 高 度 .温度 . 比 湿 和 
风速 的 预报 RMSE 分 别 小 于 20 gpm、2 K,8x10* kg， 
kg 以 及 5 ms。 

需要 说 明 的 是 ,限于 加 密探 空 实验 的 时 间 ,本 
研究 在 偏差 分 析 时 , 仅 选 择 了 2016 年 7 月 1 个 月 作 
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为 研究 时 ja] Ex 3 由 于 中 亚 地 KX F 面 植被 覆盖 相 文摘 要 集 , 2006. [Wu Xuebao, Paul Menzel, Allen Huang. Re- 
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Bias analysis and applicability evaluation of the atmospheric 
infrared sounder (AIRS) radiance in Central Asia 


MA Yufen'*, LI Ruqi’, ZHANG Meng’, Ali Mamtimin'’, ZHANG Guangxing"’ 
(1. Institute of Desert Meteorology, CMA, Urumqi 830002, Xinjiang, China; 2. Central Asia Atmospheric Science 
Research Center, Urumqi 830002, Xinjiang, China; 3. Xinjiang Meteorological Observatory, Urumqi 830002, 


Xinjiang, China) 


Abstract: Due to a scarcity of observation sites, only a small amount of conventional observation data on the 
temporal and spatial distribution of temperature and humidity in Central Asia can be obtained, which makes 
analysis difficult. High-resolution air infrared detector (AIRS) data can effectively fill the gap. In this paper, the 
radiance temperature of AIRS simulated by the input radiosonde in the Community Radiative Transfer Model was 
utilized as the reference value, deviations in the brightness temperature of the AIRS observation were analyzed, 
and the applicability of AIRS satellite data in the Central Asia numerical weather forecast operation system was 
evaluated. First, it shows that the average value of the maximum positive deviation of brightness temperature 
over the selected stations was approximately 3.3 K, and the absolute value of the maximum negative deviation 
was approximately 2.6 K. Second, the average brightness temperature of the AIRS observation in multiple 
stations was slightly higher than the overall simulated brightness temperature, and its probability density 
distribution was closer to the normal distribution curve than that of a single station. Finally, the assimilation of 
AIRS improved the prediction effect of RMAPS-CA on the geopotential height, temperature, specific humidity, 
and other high-altitude elements, but did not improve the prediction of high-altitude wind speed. For each factor, 
the assimilation improvement range of AIRS was larger at the lower and higher levels. After assimilation, the root 
mean square error of the geopotential height, temperature, specific humidity, and wind speed were less than 20 
gpm, 2 K, 8x10* kg-kg' and 5 m- s”, respectively. It should be noted that, due to the time limitations imposed on 
the encrypted sounding experiment, only one full month in July 2016 was selected as the research time period in 
the deviation analysis of this study. Due to the relatively less vegetation coverage on the underlying surface in 
Central Asia, the surface radiation was large in summer, while the area without snow cover in winter was 
relatively small. Therefore, the conclusion of the deviation analysis in this paper is not necessarily applicable to 
other seasons. In addition, considering the possibility of a business transformation of the research results, this 
study was based on Rapid-Refresh Multiscale Analysis and Prediction System-Central Asia (RMAPS-CA) to carry 
out assimilation analysis when evaluating the applicability of AIRS. The system uses the three- dimensional 
variational data assimilation (3DVAR) method. If a more advanced four-dimensional variational data assimilation 
(4DVAR) assimilation method is adopted, the assimilation effect may improve. 

Keywords: AIRS; Central Asia; radiative brightness temperature; bias analysis; applicability 


